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SUMMARY: An analysis of the relation between radio surface brightness and 
diameter, so-called £ — D relation, for planetary nebulae (PNe) is presented: i) the 
theoretical £ — D relation for the evolution of bremsstrahlung surface brightness 
is derived; ii) contrary to the results obtained earlier for the Galactic supernova 
remnant (SNR) samples, our results show that the updated sample of Galactic 
PNe does not severely suffer from volume selection effect - Malmquist bias (same 
as for the extragalactic SNR samples) and; iii) we conclude that the empirical 
£ — D relation for PNe derived in this paper is not useful for valid determination 
of distances for all observed PNe with unknown distances. 
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1. INTRODUCTION 



The relation between radio surface bright- 
nesses and diameters of supernova remnants (SNRs), 
the so-called £ — D relation, has been subject of the 
extensive discussions in the last more than fourty 
years. Due to improvements of the observational 
techniques (radio-interferometers), the several hun- 
dreds planetary nebulae (PNe) were resolved in the 
last two decades at radio frequencies, but the £ — D 
relation for PNe was not discussed until now. By 
using radio data, some statistical methods were es- 
tablished in order to determine distances to PNe. 
The main method was related to the correlation be- 
tween radius of PNe and brightness temperature - 
R — I), relation (Van de Steene and Zijlstra 1995, 
Zhang 1995, Phillips 2002). The different samples 
of Galactic PNe with known distances were defined 
in these papers. All the obtained empirical R — T& 



relations were used for determination of distances to 
PNe for which the independent distances (in order of 

R — Tf, dependence) were not obtained earlier. 

The samples of Galactic PNe are better for 
statistical analysis than the samples of Galactic 
SNRs. The selection effects should be smaller in 
the case of PN samples. However, the selection ef- 
fects surely influence the Galactic PN samples and 
the statistical determination of distances to Galactic 
PNe has to be highly uncertain. 



lowing: 



The main objectives of this paper are the fol- 

i) to derive a simple form of the theoretical 
£ — D relation for PNe by analyzing the evolution of 
radio bremsstrahlung surface brightness, 

ii) to discuss whether the updated sample of 
radio PNe is affected by the selection effects, and, 

iii) to check whether the £ — D relation is valid 
for determination of distances to PNe. 
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2. ANALYSIS AND RESULTS 

2.1. Theoretical E - D relation for PNe 

The thermal bremsstrahlung mechanism is re- 
sponsible for radiation of HII regions at radio wave- 
lengths. The bremsstrahlung volume emissivity e v of 
a PN can be shown to be (Rohlfs and Wilson 1996): 

e„[ergs s" 1 cm~ 3 Hz" 1 ] oc n 2 T~ 1/2 , (1) 

where n is the volume density and T is the thermo- 
dynamic temperature of interstellar medium (ISM). 
The surface brightness can be expressed as: 

£„ oc s v D, (2) 

where D is the diameter of PN. Combining Eqs. (1) 
and (2), we obtain: 

£„ oc n 2 T- x l 2 D. (3) 

Our next step is to express dependance of n and T 
on D. For a constant velocity mass flow the den- 
sity distribution is p = . 5 . i.e. n oc D~ x , where 
x = 2. Moreover, for the isothermal envelope with 
a power-law electron density distribution there is re- 
lationship between the shape of the density distri- 
bution and the power-law index of the radio contin- 
uum spectra (see Gruenwald and Aleman 2007, and 
references therein). Supposing that n oc D~ 2 and 
T=const. (HII regions are approximately isothermal 
at T ~ 10 4 K), we obtain the simplest form of the 
theoretical E — D relation for PNe: 

E„ oc D- 3 . (4) 

This is a standard power-law form of the E — D 
relation which can be written in general form as 
E = AD~@ , that is the same as in the case of SNRs. 

It is possible that x in density distribution is 
slightly higher, x > 2, and that the temperature is 
not strictly constant throughout the nebula. We can 
expect to see temperature gradients in PNe arising 
from radiation hardening. More energetic photons 
will travel further and when they are absorbed by the 
PN they will impart greater kinetic energy to the ions 
thereby producing a higher temperature. Using the 
numerical model results given by Evans and Dopita 
(1985), we calculate the dependence between logT e 
and logD and find the low slope (« 0.1). Therefore, 
this only slightly changes the slope of the theoretical 
E — D relation. The value (3 = 3 is then a theoretical 
lower limit, and the E — D relation could only be 
steeper, as one can see from Eq. (3). 

2.2. The empirical E — D relation for PNe 

The most important prerequisite for deriving 
a proper empirical E — D relation is defining of a 
representative sample of PNe. The distances to the 
calibrators have to be determined by accurate meth- 
ods, e.g. trigonometric or spectroscopic parallaxes of 



central stars in PNe, or by a method that uses the 
expansion of nebulae. On the other hand, all sam- 
ples suffer from the severe selection effects that arise 
from limitation in sensitivity and resolution, but the 
most severe selection effect for the Galactic samples 
of PNe is Malmquist bias; i.e. intrinsically bright 
PNe are favored because they are sampled from a 
larger spatial volume compared to any given flux lim- 
ited survey. The result is a bias against low surface 
brightness nebulae such as highly evolved old PNe. 
In this paper we use the updated sample of PNe at 
the distances less than 0.7 kpc collected by Phillips 
(2002). The influence of Malmquist bias in this sam- 
ple is limited because of the limitation in distances to 
PNe. In addition, we assume that the distances are 
accurately determined for this sample of relatively 
close PNe. The empirical E — D relation at 5 GHz 
for 44 calibrators with distances less than 0.7 kpc 
(Phillips 2002) has the form: 



E 56Hz -2.33^-10- 22 ^- 2 - 07±0 - 19 . (5) 

The parameters A and (3 are calculated by least- 
squares fitting procedure with correlation coefficient 
—0.86. The corresponding E„ — D diagram is shown 
in Fig. 1. 



1 u p 1 1 — i — i i i i 1 1 1 1 — i — i i i i 1 1 1 1 — i — i i i 1 1. 




-| q 25 | i i i I i i i I i i i 

0.01 0.1 1 10 

D[pc] 

Fig. 1. The E — D diagram at 5 GHz for 44 Galac- 
tic PNe with distances less than 0.7 kpc. 

The form of Eq. (5) is very close to the so- 
called trivial E — D form with (3 = 2 (for details see 
Arbutina et al. 2004). The additional test in order to 
estimate the validity of Eq. (5) pertains to the possi- 
ble dependence between the luminosity and diameter 
of PNe. The L u — D diagram is shown in Fig. 2. The 
scatter in L v — D plane shows that the correlation 
between L v and D is poor (correlation coefficient = - 
0.06) and therefore the physical dependence between 
L and D could not be confirmed by this statistical 
procedure. 
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Fig. 2. The L- D plot at 5 GHz for U Galactic 
PNe with distances less than 0. 7 kpc. 



3. DISCUSSION 



The theoretical £„ — D relation (Eq. (4)) for 
PNe, derived in this paper, describes a trend of de- 
creasing radio surface brightness with increasing di- 
ameter of an object. The radiation mechanism used 
in this simple derivation is thermal bremsstrahhmg. 
This is the basic process of production of the radio 
radiation in HII regions. The theoretically derived 
slope (/? = 3) is steeper than the slope from the em- 
pirical relation given by Eq. (5). This discrepancy 
can be explained by the low quality of the sample of 
Galactic PNe or by the assumptions used in deriva- 
tion of theoretical relation. Due to small variation 
in power-law density distribution with x > 2 (Gru- 
enwald and Aleman 2007, and references therein) 
and approximately constant temperature of expand- 
ing envelope of PNe, theoretical slope can be slightly 
steeper than in Eq. (4). Therefore, we conclude 
that the theoretical relation has the correct form, 
but our empirical relation is under influence of bi- 
ases that could make the slope shallower. On the 
other hand, there are some attempts to show that 
evolution of PNe are not linear in log-log scales (e.g. 
Phillips 2004). These different dependences cannot 
be derived from the thermal bremsstrahlung radia- 
tion formula (Eq. (1)). 

A very interesting feature regarding the em- 
pirical relation for Galactic PNe (Eq. (5)) is that 
the slope is approximately equal to the slope of triv- 
ial X — D relation. Therefore, we conclude that 
Malmquist bias is not so severe as in cases of Galac- 



tic SNR samples. This slope (j3 « 2) was obtained 
for the extragalactic samples of SNRs (except M82 
sample) where Malmquist bias is small, because all 
the SNRs are at the approximately same distance 
(see Urosevic 2002, Urosevic et al. 2005). 

The large scatter in L v — D plane (Fig. 2) 
suggests that the slope in Eq. (5) docs not have 
real and valid physical interpretation. It is a kind 
of luminosity-diameter scattering artefact which pro- 
duces the trivial X oc D~ 2 form. Therefore, the rela- 
tion defined by Eq. (5) is not precise enough for de- 
termination of valid distances to Galactic PNe. This 
is due to the different biases: the limitations in sen- 
sitivity and resolution of radio surveys, the source 
confusion, Malmquist bias (in mild form), mixture 
of different types of PNe in the same sample, and 
insufficient precision in determining the distances to 
the 44 calibrators. 



4. SUMMARY 



The main results of this paper may be sum- 
marized as follows: 

i) The theoretical X„ — D relation for the radio 
evolution of thermal bremsstrahlung surface 
brightness of PNe in form of X„ oc D~ 3 is de- 
rived. 

ii) Our results show that the updated sample of 
Galactic PNe does not severely suffer from vol- 
ume selection effect - Malmquist bias (same 
as in cases of the extragalactic SNR samples) . 
This is opposite to results obtained earlier for 
the Galactic SNR samples. 

iii) Due to analysis of the L u — D dependence, we 
conclude that the X„ — D relation for Galactic 
PNe is not useful for reliable determination of 
distances for all observed PNe with unknown 
distances. 

The above observation leads to the more gen- 
eral comment that PNe may have very different ini- 
tial conditions leading to independent evolutionary 
paths. These paths could follow the same theoreti- 
cal X — D curve but with varying intercepts, leading 
to the scatter such as the one found in this paper. 
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IlpHKa3aHa je aHajnoa t3b. X — D pc- 
Jiairnje H3Mefjy noBpniHHCKor cjaja Ha pa^HO- 
(JjpeKBeHiTHjaMa h ^.njaMCTpa njiaHCTapHHx 
MarjiHHa (IIM): i) vL3BeAen& je TeopnjcKa £ — D 
pejiairnja 3a eBOJiyirnjy noBpniHHCKor cjaja 
CTBopeHor 3aKOMHHM 3pa^ieH)eM; ii) cynpoTHO 
pe3yjiTaTHMa ^oSnjeHHM pamije 3a y3opKC 
ca^HBbCHe op, TajiaKTEraKHx ociaiaKa cyncp- 
hobhx, HaniH pe3yjiTaTH noKa3yjy aa HajHOBujc 



(|jopMHpaHH y3opaK TajiaKTH^KHx IIM He Tpnn 
bcjihkh yTHn,aj 36or 3anpcMHHCKor cejiCKiiHOHor 
e(})eKTa, T3B. MajiMKBHCTOBor cejiCKHHOHor 
e(})eKTa (hcto Baacn 3a BaHrajiaKTiraKe y3opaKe 
ocTaTaKa cynepHOBHx); h iii) 3aKJi>y*ryjeMO p.& 
CMnnpHjcKa £ — - D pejiairnja 3a IIM n3Be,neHa 
y obom pa,ay Hnje ynoTpe6jLHBa 3a noy3^aHa 
o.nper)iiBaH>a j\a,jbvm& p,o cbhx nocMaTpamix IIM 
ca Heno3HaTHM .zjajLHHaMa. 
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